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ABSTRACT:  The  majority  of  cephalopods  are  thought  to  have  limitations  arising  from  physiology  and 
locomotion  that  exclude  them  from  shallow,  highly  variable,  euryhaline  environments.  The  brief 
squid  Lolliguncula  brevis  may  be  a  notable  exception  because  it  tolerates  low  salinities,  withstands  a 
wide  range  of  environmental  conditions,  and  swims  readily  in  shallow  water.  Little  is  known  about 
the  distribution  of  L.  brevis  in  Chesapeake  Bay,  a  diverse  and  highly  variable  estuary.  Therefore,  a 
survey  of  L.  brevis  was  conducted  in  the  Virginia  portion  of  Chesapeake  Bay  from  1993  to  1997  using 
a  9.1  m  otter  trawl,  and  the  effects  of  selected  factors  on  squid  presence  were  assessed  using  logistic 
regression  analysis.  During  spring  through  fall,  L.  brevis  was  collected  over  a  wide  range  of  bottom- 
water  salinities  (17.9  to  35.0%o),  bottom-water  temperatures  (8.1  to  29.6°C),  bottom-water  dissolved 
oxygen  levels  (1.9  to  14.6  mg  O2  1“*),  and  depths  (1.8  to  29.9  m),  but  it  was  not  present  in  trawls  con¬ 
ducted  during  winter  in.  L.  brevis,  especially  juveniles  <  60  mm  dorsal  mantle  length  (DML),  were 
abundant,  frequently  ranking  in  the  upper  12  %  of  overall  annual  nektonic  trawl  catches,  and  during 
the  fall  of  some  years,  ranking  second  to  anchovies.  The  probability  of  catching  a  squid  increased  in 
Chesapeake  Bay  at  higher  salinities  and  water  temperatures,  and  was  much  greater  in  normoxic  than 
in  hypoxic  waters;  these  variables  had  a  profound  influence  on  both  annual  and  seasonal  variability 
in  distribution.  Salinity  had  the  largest  influence  on  squid  distribution,  with  squid  being  completely 
absent  from  the  bay  when  salinity  was  <17.9%o  and  most  abundant  in  the  fall  when  salinity  was  high¬ 
est  (despite  declines  in  water  temperature).  Squid  were  most  prevalent  at  depths  between  10  and 
15  m.  The  results  of  this  study  suggest  that  L.  brevis  is  an  important  component  of  the  Chesapeake 
Bay  ecosystem  when  salinities  and  water  temperatures  are  within  tolerance  limits  and  that  unlike 
other  squids,  L.  brevis  may  be  well-equipped  for  an  inshore,  euryhaline  existence. 
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INTRODUCTION 

Most  squids  are  unable  to  osmoregulate  and  tolerate 
low  salinities  (Hendrix  et  al.  1981,  Mangum  1991),  pre¬ 
venting  the  vast  majority  of  squids  from  invading  fresh 
or  estuarine  waters  —  2  significant  aquatic  habitats 
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where  cephalopods  are  poorly  represented.  One 
notable  squid,  the  brief  squid  Lolliguncula  brevis,  is 
the  only  species  of  cephalopod  frequently  found  in 
low-salinity  estuaries  (Vecchione  1991a),  where  it 
tolerates  salinities  as  low  as  8.5 %o  for  brief  periods 
(Laughlin  &  Livingston  1982).  Within  shallow  bays  and 
estuaries,  there  is  evidence  that  L.  brevis  withstands 
low  dissolved  oxygen  levels  (Vecchione  1991b)  and  a 
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wider  range  of  water  temperatures  than  most  cephalo- 
pods  (Hixon  1980,  Roper  et  al.  1984).  This  species  is 
smaller  than  most  loliginids,  seldom  reaches  sizes 
greater  than  110  mm  in  dorsal  mantle  length  (DML) 
(Hixon  1980,  Vecchione  et  al.  1989),  has  rounded  wide 
fins,  and  appears  to  swim  at  low  velocities  (Finke  et  al. 
1996,  Bartol  et  al.  2001a,  b).  These  traits  are  beneficial 
for  maneuvering  in  nearshore  waters  of  the  Gulf  of 
Mexico  and  along  the  western  Atlantic  coast  from 
Delaware  to  Brazil,  where  L.  brevis  is  commonly  found 
(Vecchione  et  al.  1989). 

Little  is  known  about  Lolliguncula  brevis  distribution 
in  Chesapeake  Bay.  This  is  surprising,  given  that  L. 
brevis  presumably  has  exceptional  tolerance  among 
cephalopods  to  salinity  and  other  physical  conditions 
and  that  Chesapeake  Bay  is  one  of  the  largest,  most 
diverse,  physically  variable  estuaries  in  the  world,  with 
tributaries  and  rivers  within  the  system  constituting  a 
shoreline  a  few  thousand  kilometers  long.  Within 
Chesapeake  Bay,  physical  and  biological  processes 
operate  rapidly  and  intensely  on  many  different  tem¬ 
poral  and  spatial  scales:  large-scale  pycnoclines 
extending  almost  the  entire  length  of  the  bay  fre¬ 
quently  develop  throughout  the  spring  (Brandt  et  al. 
1986);  significant  phytoplankton  blooms  arise  in  the 
late  spring/early  summer  (Pinet  1992);  regions  of  bot¬ 
tom-water  hypoxia  may  develop  in  the  late  summer  for 
weeks  or  even  between  tides  when  organic  matter 
deposition  is  high  and  vertical  mixing  is  difficult 
because  of  water-column  stratification  (Webb  &  D'Elia 
1980);  hundreds  of  northern  and  southern  species  of 
fishes  converge  in  the  bay  throughout  the  year  (Geer  & 
Austin  1995);  and  dramatic  decreases  in  salinity  and 
water  temperature  occur  rapidly  in  the  spring  and  win¬ 
ter,  respectively  (Brandt  et  al.  1986).  Given  that  most 
squids  evolved  to  reside  in  physically  stable,  offshore 
environments  where  interactions  with  fishes  are  low 
relative  to  inshore  environments  (O'Dor  &  Webber 
1986,  Wells  1994),  the  presence  of  L.  brevis  in  a  highly 
variable,  euryhaline  environment  rich  in  nektonic 
fauna  such  as  Chesapeake  Bay  is  intriguing.  Under¬ 
standing  the  distribution  of  L.  brevis  in  Chesapeake 
Bay  in  relation  to  physical  factors  promises  insight  into 
the  ecology  of  a  unique  squid  that  lives  in  an  environ¬ 
ment  that  tests  its  physiological  limitations. 

To  learn  about  the  distribution  of  brief  squid  in 
Chesapeake  Bay,  a  survey  was  conducted  in  the  Vir¬ 
ginia  portion  of  the  bay  from  1993  to  1997  using  an 
otter  trawl.  The  primary  objective  of  this  study  was  to 
determine  what  and  how  selected  abiotic  factors  affect 
brief  squid  presence  within  Chesapeake  Bay.  The  abi¬ 
otic  factors  of  greatest  interest  were:  tidal  stage,  bot¬ 
tom  depth,  salinity,  water  temperature,  dissolved  oxy¬ 
gen  and  water  clarity.  Since  temporal  variation  in  the 
Chesapeake  Bay  ecosystem  is  high,  the  effects  of  year 


and  season  also  were  considered.  Two  secondary  ob¬ 
jectives  of  the  study  were  to  provide  some  assessment 
of  the  abundance  of  brief  squid  in  Chesapeake  Bay  rel¬ 
ative  to  other  nekton  and  to  document  the  size  distrib¬ 
ution  of  brief  squid  in  Chesapeake  Bay  throughout  the 
year. 

MATERIALS  AND  METHODS 

Sampling  design.  Data  on  Lolliguncula  brevis  distri¬ 
bution  were  collected  from  1993  to  1997  during  the 
Virginia  Institute  of  Marine  Science  (VIMS)  Trawl  Sur¬ 
vey  Monitoring  Program,  which  was  established  to 
generate  annual  indices  for  juvenile  marine  and  estu¬ 
arine  finfish  and  invertebrates  within  Virginia  waters. 
During  this  study,  both  the  mainstem  of  the  lower 
Chesapeake  Bay  from  the  Bay  mouth  to  37° 40'  N  and  3 
major  Virginia  tributaries,  the  James,  York,  and  Rap- 
pahanock  Rivers,  were  sampled  (see  Fig.  1).  The  gear 
used  consisted  of  a  9.1  m  semi-balloon  otter  trawl  with 
3.8  cm  body  mesh,  3.2  cm  cod-end  mesh,  1.3  cm  cod- 
end  liner  mesh,  a  tickler  chain,  and  steel  China-V 
doors  (71  cm  x  48  cm).  Each  tow  lasted  5  min  (bottom 
time)  and  was  performed  at  a  speed  of  1.29  m  s“^. 

A  stratified  random  design  with  stratification  by  lati¬ 
tude  and  water  depth  was  used  to  sample  sites  within 
the  Chesapeake  Bay  mainstem.  Three  latitudinal  strata 
(36°55'  -  37°10'N;  37°10'  -  37°25'N;  and  37°25'  - 
37°40'N)  and  4  depth  strata  (eastern  bay  shoal  areas 
[3.6-9. 1  mj;  western  bay  shoal  areas  [3.6-9. 1  mj;  plain 
areas  [9.1-13  mj;  and  deep  areas  [>13  mj)  were  consid¬ 
ered,  resulting  in  a  total  of  12  strata.  The  strata  were 
fixed  in  time  and  did  not  change  with  tidal  cycles.  Sam¬ 
pling  sites  were  selected  randomly  within  the  12  strata 
using  computerized  files  of  the  National  Ocean  Sur¬ 
vey's  Chesapeake  Bay  bathymetric  grid  system.  Sam¬ 
pling  within  the  Chesapeake  Bay  mainstem  was  con¬ 
ducted  monthly  from  April  through  December  and  in 
February  of  each  year.  Within  each  of  the  3  latitudinal 
strata  for  every  month  considered  during  a  given  year, 
sampling  included  4  stations  within  the  plain  stratum,  3 
stations  within  the  deep  stratum,  and  either  3  (in  May 
through  November)  or  2  (in  December,  February,  and 
April)  stations  within  each  of  the  2  shoal  strata.  There¬ 
fore,  33  to  39  stations  were  sampled  each  month  (except 
January  and  March  when  no  samples  were  collected) 
within  the  mainstem  of  Chesapeake  Bay. 

Within  the  3  tributaries,  fixed  stations  located  in  the 
center  of  the  river  channels  at  8  km  intervals  from  the 
river  mouths  to  the  freshwater  interfaces  were  sam¬ 
pled  monthly  from  1993  to  1997.  In  1996  and  1997, 
additional  stations  were  sampled  using  a  stratified  ran¬ 
dom  design  partitioned  by  region  and  depth.  Each 
river  was  divided  into  4  regions,  but  the  number  of 
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depth  strata  (2  to  5)  within  each  region  varied  accord¬ 
ing  to  the  dimensions  of  the  tributaries.  Generally,  25 
to  62  tributary  stations  were  sampled  each  month 
within  the  3  tributaries. 

Bottom  depth,  bottom-water  salinity,  bottom-water 
temperature,  bottom- water  dissolved  oxygen,  and  Sec- 
chi  depth  (a  measure  of  water  clarity)  were  recorded  at 
all  stations  immediately  following  each  tow.  Bottom 
depth  (±0.1  m)  was  measured  using  an  ultrasonic 
depth  sounder,  salinity  (±0.1  %o),  water  temperature 
(±0.1°C)  and  dissolved  oxygen  (±0.1  mg  O2  1“*)  were 
measured  using  a  Hydrolab  multiprobe  (Hydrolab 
Corporation,  Austin,  Texas);  and  Secchi  depth  (±0.1  m) 
was  measured  using  a  Secchi  disk.  Tidal  stage  was 
estimated  using  Chesapeake  tidewater  tide  logs 
(Pacific  Publishers,  Bolinas,  California)  and  Tide  and 
Currents  Version  2.0  software  (Nautical  Software,  Inc., 
Beaverton,  OR).  Latitude  and  longitude  at  the  begin¬ 
ning  and  ending  of  each  tow  were  determined  using 
Loran  C  conversions  or  a  Global  Positioning  System 
(GPS)  receiver.  Tow  distance  was  defined  as  the  dis¬ 
tance  between  beginning  and  ending  latitude/longi¬ 
tude  coordinates. 

The  number  and  DMLs  of  squid  and  the  number  and 
types  of  vertebrate  nekton  captured  in  the  trawls  were 
recorded.  Occasionally,  squid  were  dissected  to  char¬ 
acterize  the  approximate  stage  of  gonad  maturity. 
Male  squid  were  considered  mature  when  they  had  a 
well-developed  penis,  tightly  packed  spermatophores 
within  the  spermatophoric  sac,  and  a  yellowish  or 
whitish  testis.  Female  squid  were  considered  ripe 
when  their  internal  oviduct  was  filled  with  mature, 
amber-colored  eggs. 

One  limitation  of  the  trawl  survey  data  was  that 
squid  captured  prior  to  1997  in  Chesapeake  Bay  were 
simply  classified  as  squid  and  assumed  to  be  LoIIigun¬ 
cula  brevis.  Examinations  of  squid  captured  in  1997 
and  1998  revealed  that  several  large  (>110  mm  DML) 
squid  captured  in  trawls  at  the  Chesapeake  Bay  mouth 
in  high  salinities  were  Loligo  pealei,  a  loliginid  that  is 
typically  larger  than  LoIIiguncula  brevis,  reaching 
sizes  up  to  500  mm  DML  (Vecchione  et  al.  1989).  The 
number  of  squid  identified  in  1997  and  1998  as  Loligo 
pealei  was  small  relative  to  those  identified  as  LoIIi¬ 
guncula  brevis,  however,  and  no  squid  less  than  110 
mm  has  been  identified  to  date  as  Loligo  pealei.  Since 
LoIIiguncula  brevis  rarely  grows  larger  than  110  mm 
DML  (Hixon  1980),  all  squid  >110  mm  DML  were  con¬ 
sidered  to  be  Loligo  pealei  and  eliminated  from  the 
analysis.  Conversely,  all  squid  <110  mm  were  consid¬ 
ered  to  be  LoIIiguncula  brevis  and  included  in  the 
analysis.  Although  this  does  not  guarantee  that  only  L. 
brevis  were  considered  in  this  study,  the  number  of 
Loligo  pealei  in  the  samples  was  probably  negligible. 
This  assertion  is  based  on:  (1)  independent  trawling 


performed  in  the  Chesapeake  Bay  ecosystem  for  pur¬ 
poses  of  collecting  specimens  for  separate  swimming 
physiology  studies,  where  1  out  of  1243  squid  captured 
were  L.  pealei-,  (2)  identifications  of  present  and  past 
trawl-survey  specimens,  which  reveal  small  numbers 
of  L.  pealei,  and  (3)  the  physiology  and  ecology  of  L. 
pealei  (L.  pealei  avoids  euryhaline  environments,  pre¬ 
ferring  areas  >32%o  and  >40  m  [Hixon  1980],  which 
are  rare  in  Chesapeake  Bay). 

Statistical  analysis  to  assess  the  effects  of  selected 
factors  on  squid  presence.  Since  the  primary  objective 
of  this  study  was  to  determine  what  and  how  selected 
physical  variables  influence  squid  presence  within 
Chesapeake  Bay,  the  outcome  or  dependent  variable 
was  simply  the  presence  or  absence  of  squid  at  each 
sampling  site  and  the  independent  variables  were  the 
selected  factors  of  interest  (e.g.  depth,  salinity,  etc.). 
Logistic  regression  using  the  logit  link,  a  standard 
method  of  analyzing  dichotomous,  discrete  outcome 
variables  (Hosmer  &  Lemeshow  1989),  was  used  to 
model  the  presence/absence  of  brief  squid  in  trawl 
catches.  The  specific  form  of  the  logistic  regression 
model  is  as  follows: 

_  exp((x  +  PiXi  +  ^2X2 . . 

~  l  +  exp(a  +  [3iXi+p2X2...p,Xi) 

where  p  =  event  probability  (i.e.,  probability  of  catch¬ 
ing  at  least  1  squid),  a  is  a  constant,  X,  is  an  indepen¬ 
dent  variable,  i  (e.g.  salinity,  temperature,  etc.),  and  Pj 
is  a  coefficient  of  the  independent  variable,  Xj. 

The  coefficients  (P)  of  the  variables  included  in  the 
model  are  most  interpretable  when  expressed  in  terms 
of  the  odds  ratio  (\|/):  \|/  =  e^.  Both  scaled  continuous  and 
discrete  variables  are  considered  in  the  logistic  regres¬ 
sion  model.  For  scaled  continuous  independent  vari¬ 
ables,  the  odds  ratio  (y)  is  a  measure  of  how  much 
more  likely  (or  unlikely  if  the  odds  ratio  is  <1)  it  is  for 
the  outcome  (in  this  case,  the  presence  of  at  least  1 
squid  in  a  trawl)  to  be  present  when  the  variable  is 
increased  by  1  unit.  Consequently,  the  interpretation 
of  a  continuous  independent  variable  within  a  logistic 
regression  model  will  depend  heavily  on  what  scaling 
units  (e.g.  2%o  or  2°C)  are  selected  for  the  variable,  and 
thus  selected  units  should  be  of  biological  relevance 
(Hosmer  &  Lemeshow  1989).  For  example,  it  may  be 
more  biologically  informative  to  determine  how  catch 
probabilities  vary  with  every  2°C  change  in  water  tem¬ 
perature  rather  than  with  every  0.2°C  change  when 
ecological  variations  may  be  very  slight.  Therefore, 
continuous  variables  were  scaled  and  coded  according 
to  units  of  reasonable  biological  relevance  and  to 
ensure  there  was  a  sufficient  sample  size  in  each  cate¬ 
gory  level.  Continuous  variables  considered  in  this 
analysis  with  their  respective  scaling  units  denoted  in 
parentheses  were:  bottom  depth  (2  m),  bottom  salinity 
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(2%o),  bottom  temperature  (2°C),  bottom  dissolved  oxy¬ 
gen  (3  mg  O2  r^),  tow  distance  (100  m)  and  Secchi 
depth  (0.5  m). 

For  discrete  independent  variables,  the  odds  ratio  is 
a  measure  of  how  much  more  likely  (or  unlikely)  it  is 
for  the  outcome  to  be  present  relative  to  a  reference 
category  (Hosmer  &  Lemeshow  1989,  Agresti  1990). 
For  example,  within  a  discrete  variable  such  as  season, 
the  coefficients  of  fall,  summer,  and  winter  may  be 
expressed  relative  to  the  coefficient  of  spring.  The  dis¬ 
crete  variables  considered  in  this  analysis  were:  year 
(1993,  1994,  1995,  1996,  and  1997),  season  (spring,  fall, 
summer,  and  winter),  and  tidal  stage  (early  flood,  max¬ 
imum  flood,  late  flood,  slack  before  ebb,  early  ebb, 
maximum  ebb,  late  ebb,  and  slack  before  flood). 

Two  significant  findings  were  detected  during  the 
early  stages  of  logistic  regression  model  construction: 

(1)  no  squid  were  caught  in  trawls  during  winter,  and 

(2)  no  squid  were  caught  in  areas  with  salinities  below 
17.5%o.  Since  there  was  a  zero  probability  of  catching 
squid  in  areas  sampled  during  the  winter  or  in  waters 
<17.5%o,  all  sites  meeting  these  criteria  were  elimi¬ 
nated  and  not  used  in  logistic  regression  models.  A 
total  of  2762  sampling  sites  were  eliminated,  leaving 
2393  sampling  sites  for  consideration  within  the  logis¬ 
tic  regression  models. 

Logistic  regression  model  construction  and  refine¬ 
ment  were  performed  following  procedures  suggested 
by  Hosmer  &  Lemeshow  (1989:  Chapter  4).  First,  uni¬ 
variate  logistic  regression  analyses  were  performed  on 
each  variable.  Significance  was  assessed  at  a  =  0.25 
using  the  -2  log  likelihood  estimation  chi-square,  and 
non-significant  variables  were  eliminated  from  consid¬ 
eration.  This  level  of  significance,  rather  than  the  con¬ 
ventional  a  =  0.05  level,  was  selected  to  minimize  the 
chance  of  excluding  variables  that  might  be  important 
outcome  predictors  when  considered  collectively  with 
other  variables  in  a  multivariate  analysis  (Mickey  & 
Greenland  1989).  At  the  univariate  stage,  scaling  of 
discrete  variables  also  was  examined;  when  several 
categories  within  a  discrete  variable  had  similar  coeffi¬ 
cients  and  odds  ratios,  they  were  combined  to  form 
broader  categories  (Hosmer  &  Lemeshow  1989). 

Upon  completion  of  the  univariate  analyses,  all  sig¬ 
nificant  variables  were  included  in  a  multivariate 
model.  Variables  that  were  not  significant  at  the  a  = 
0.05  level  within  the  multivariate  model  and  that  did 
not  contribute  to  wide  fluctuations  in  coefficients  after 
removal  were  eliminated,  and  further  multivariate 
logistic  regression  models  were  performed.  Once  all 
essential  variables  were  determined,  the  scale  of  con¬ 
tinuous  variables  was  examined  so  that  assumptions  of 
linearity  in  the  logit  were  not  violated.  Linearity  of 
scaled  continuous  variables  was  examined  by  dividing 
values  into  quartiles  (or  smaller  groupings  when 


higher  resolution  was  necessary)  and  treating  continu¬ 
ous  variables  as  discrete  variables  in  the  multivariate 
model  using  the  lowest  quartile  (or  group)  as  the  refer¬ 
ence  level.  This  procedure  revealed  that  it  was  some¬ 
times  necessary  to  truncate  the  range  of  categories 
within  a  variable,  increase  unit  size,  and  treat  some 
scaled  continuous  variables  as  discrete  and  nominal 
when  assumptions  of  linearity  were  violated.  Although 
this  scaling  may  reduce  resolution  within  the  model,  it 
stabilizes  parameter  estimates,  produces  goodness-of- 
fit  measures  of  greater  reliability,  and  makes  the  model 
more  interpretable  (Hosmer  &  Lemeshow  1989, 
Agresti  1990). 

Once  all  relevant  variables  were  determined  and 
expressed  in  the  correct  scale,  further  model  refine¬ 
ment  involved  the  consideration  of  interactions. 
Because  of  practical  considerations,  it  is  not  possible  to 
include  all  interactions  in  logistic  regression  models 
that  have  many  variables.  Only  those  interactions  that 
make  biological  sense  should  be  investigated  (Hosmer 
&  Lemeshow  1989).  Thus,  we  restricted  our  search  for 
interactions  to  those  involving  the  4  abiotic  factors  of 
greatest  interest:  depth,  salinity,  water  temperature, 
and  dissolved  oxygen  level.  Each  interaction  involving 
these  factors  was  added  separately  to  the  final,  main- 
effects-only,  multivariate  model.  A  determination  as  to 
whether  to  include  the  interaction  in  further  models 
was  based  on  both  the  significance  of  the  -2  log  likeli¬ 
hood  estimation  of  the  new  model  and  the  p-value  of 
the  interaction.  Interactions  that  were  significant  at  a  = 
0.05  and  that  provided  significant  improvement  over 
the  main-effects-only  model  were  added  to  the  final 
model.  To  decouple  the  effects  of  factors  involved  in 
significant  interactions,  logits  for  all  factor  combina¬ 
tions  were  computed  and  significance  assessed  at 
a  =  0.05  using  procedures  described  in  detail  on 
p.  102-106  of  Hosmer  &  Lemeshow.  The  fit  of  the  final, 
most  parsimonious  model  was  assessed  using  Hosmer 
and  Lemeshow  goodness-of-fit  tests  and  various  logis¬ 
tic  regression  diagnostics  (Hosmer  &  Lemeshow  1989: 
p  149-170).  There  was  no  overdispersion  in  the  model, 
i.e.,  data  did  not  deviate  significantly  from  the  ex¬ 
pected  binomial  distribution. 

NMFS/NEFSC  bottom  trawl  surveys.  Size  distri¬ 
bution  and  catch  location  data  on  Lolliguncula  brevis 
collected  during  spring  and  autumn  bottom-trawl 
surveys  by  the  National  Marine  Fisheries  Service 
(NMFS)/Northeast  Fisheries  Science  Center  (NEFSC) 
for  1993  to  1997  were  considered  in  this  study.  These 
surveys  were  conducted  within  coastal  waters  (30  to 
1200  m)  of  the  Atlantic  Ocean  from  Cape  Hatteras, 
North  Carolina,  to  the  Gulf  of  Maine.  Approximately 
300,  30  min  tows  were  performed  biannually  at  ran¬ 
domly  selected  sites  using  a  #36  otter  trawl  with  a 
1.26  cm  cod-end  mesh  liner.  Dorsal  mantle  lengths  of 
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L.  brevis  captured  in  each  tow  were  recorded,  and 
capture  location  and  depth  were  noted.  Further  infor¬ 
mation  on  NMFS/NEFSC  bottom-trawl  procedures 
may  be  found  in  'Fisherman's  Reports'  for  1993 
to  1997  (NMFS,  NEFSC,  Woods  Flole,  Massachusetts 
02543)  and  at:  http://www.nefsc.nmfs.gov/sos/vesurv 
/vesurv.html. 


RESULTS 

LoIIiguncula  brevis  was  captured  at 
545  of  the  5155  sampling  sites.  In 
trawls  containing  squid,  the  mode  was 
1  squid  caught  per  trawl,  the  median 
was  4  squid  caught  per  trawl,  the 
mean  was  16.2  ±  31.9  (SD)  squid 
caught  per  trawl,  and  the  range  was  1 
to  325  squid  caught  per  trawl.  The 
majority  of  L.  brevis  captured  during 
the  survey  were  collected  in  the  main- 
stem  of  Chesapeake  Bay  and  not  the 
tributaries,  and  the  greatest  overall 
number  of  squid  collected  in  trawls 
occurred  within  the  central  channel  of 
the  lower  bay  and  along  the  eastern 
portion  of  the  bay  (Eig.  1).  L.  brevis 
were  captured  during  spring,  summer, 
and  fall,  but  no  squid  were  captured  in 
trawls  conducted  during  winter.  Estu¬ 
arine  waters  of  0.9  to  43.9  m,  0.0  to 
37.1%o,  0.3  to  31.0°C,  0.1  to  14.6  mg 
O2  r^,  and  were  sampled;  brief  squid 
were  collected  over  a  wide  range  of 
depths  (1.8  to  29.9  m),  bottom-water 
salinities  (17.9  to  35.0  %o),  bottom- 
water  temperatures  (8.1  to  29.6°C), 
and  bottom-water  dissolved  oxygen 
levels  (1.9  to  14.6  mg  O2  1”^). 

In  terms  of  overall  numbers  of  ani¬ 
mals  caught  annually  in  trawls,  LoIIi¬ 
guncula  brevis  ranked  9  out  of  94  total 
fish/squid  species  in  1993  (90th  per¬ 
centile),  18  out  of  92  total  fish/squid 
species  in  1994  (80th  percentile),  6  out 
of  94  total  fish/squid  species  in  1995 
(94th  percentile),  46  out  of  98  total 
fish/squid  species  in  1996  (53rd  per¬ 
centile),  and  12  out  of  103  total  fish/ 
squid  species  in  1997  (88th  percentile). 
When  distinctions  are  made  between 
Chesapeake  Bay  mainstem  and  tribu¬ 
tary  sites,  it  is  clear  that  squid  ranked 
consistently  higher  annually  within 
the  Chesapeake  Bay  mainstem  than 


within  the  tributaries  (Eig.  2).  Moreover,  squid  gener¬ 
ally  ranked  higher  in  the  summer  and  fall  than  in  the 
winter  and  spring  both  in  the  Chesapeake  Bay  main- 
stem  and  the  tributaries.  Brief  squid  were  especially 
abundant  within  the  Chesapeake  Bay  mainstem  dur¬ 
ing  the  fall  of  1995,  when  squid  ranked  second  (behind 
bay  anchovies  Anchoa  mitchilli]  out  of  53  (96th  per¬ 
centile)  squid/fish  caught  (Eig.  2). 

Most  of  the  squid  captured  during  this  survey  were 
<60  mm  DML,  although  some  larger  squid  (>70  mm 
DML)  were  present  in  the  spring  of  most  years  (i.e., 


Fig.  1.  Map  of  Chesapeake  Bay  study  site.  Mainstem  of  lower  Chesapeake  Bay 
from  the  bay  mouth  (southern  sampling  limit)  to  37°40'N  (northern  sampling 
limit)  and  3  Virginia  tributaries,  the  James,  York,  and  Rappahanock  Rivers, 
were  sampled.  Upstream  and  downstream  tributary  boundaries  are  depicted. 
Number  of  brief  squid  LoIIiguncula  brevis  captured  in  trawls  from  1993  to  1997 
in  Chesapeake  Bay  are  shown  (circles)  as  well  as  locations  where  egg  capsules 
were  found  (x).  All  squid  in  this  figure  were  <110  mm  in  dorsal  mantle  length. 
Each  tow  lasted  5  min  (bottom  time)  and  was  performed  at  a  speed  of  9.7  km^'h 
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1993,  1995,  1996,  1997)  and  comprised  a  notice¬ 
able  proportion  of  spring  catches  (Fig.  3).  The 
mean  DML  of  squid  decreased  from  spring  to 
summer  in  each  year  of  the  study.  However, 
with  the  exception  of  1994,  when  size  increased 
substantially  from  summer  to  fall,  there  was  lit¬ 
tle  difference  in  the  mean  DML  of  squid  caught 
in  summer  and  fall.  In  contrast,  a  clear  increase 
in  Lolliguncula  brevis  mantle  length  from 
spring  to  fall  was  observed  in  coastal  waters 
from  Delaware  to  Florida  in  NMFS/NEFSC  sur¬ 
veys.  In  these  surveys,  mean  L.  brevis  DML  was 
32  mm  DML  (±12  SD)  during  the  spring  and  62 
mm  DML  (±23  SD)  during  the  fall  (Fig.  4).  In 
the  NMFS/NEFSC  surveys,  L.  brevis  were 
caught  in  depths  shallower  than  36  m,  with  the 


Fig.  3.  Lolliguncula  brevis.  Length/frequency  data  of 
brief  squid  captured  in  spring,  summer,  and  fall  dur¬ 
ing  fhe  VIMS  trawl  surveys  from  1993  to  1997. 
Annual  data  as  well  as  cumulative  data  for  the  entire 
sampling  period  (lower  right  graph)  are  depicted.  For 
each  of  the  3  seasons,  proportion  of  total  catches 
within  each  of  11  size  classes  is  presented.  Total  num¬ 
ber  (n)  and  mean  dorsal  length  (m)  of  squid  caughf 
within  the  3  seasons  are  included  above  the  graphs 


O 

CO 


Fig.  2.  Lolliguncula  brevis.  Seasonal  (winter,  spring,  sum¬ 
mer,  and  fall)  and  annual  (year)  rankings  for  brief  squid 
caught  by  trawl  in  mainstem  and  tributary  sections  of 
Chesapeake  Bay  in  1993  to  1997.  Rankings  are  expressed  as 
percentiles  and  are  based  on  the  total  number  of  brief  squid 
caught  relative  to  the  total  number  of  each  fish  species 
caught  per  season  or  per  year  at  mainstem  and  tributary 
regions  of  Chesapeake  Bay.  Fractions  above  bars  denote  the 
exact  rank  of  squid  out  of  the  total  number  of  squid/fish  spe¬ 
cies  considered:  lower  fractions  are  results  from  mainstem 
stations,  upper  fractions  for  tribufary  sfations 


majority  of  squid  being  caught  in  waters  shallower 
than  23  m. 

Tow  distance  and  Secchi  depth  did  not  significantly 
influence  the  probability  of  catching  at  least  1  squid 
in  a  trawl,  and  were  thus  eliminated  from  the  final 
logistic  regression  model.  The  factors  and  respective 
categories  considered  in  the  final,  most  parsimonious 
model  are  shown  in  Table  1.  The  logistic  regression 
model  was  highly  significant  (-2  log-likelihood  chi 
square  =-786.15;  goodness-of-fit  statistic,  G  =  665.28; 
df  =  15;  p  <0.001),  and  the  model  adequately  fit  the 
data  (Hosmer  &  Lemeshow  test,  p  =  0.359)  (Table  1). 
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Fig.  4.  LoIIiguncula  brevis.  Length/frequency  data  of  squid 
collected  during  spring  and  fall  NMFS/NEFSC  surveys  (1993 
to  1997)  along  the  coastline  of  the  US  from  Delaware  to 
Florida.  Proportion  of  total  spring  and  fall  catches  within  each 
of  10  size  classes  is  presented.  In  the  fall  survey  518  squid 
were  caught,  in  the  spring  survey  713  squid 

Tidal  stage  and  depth  affected  capture  probability. 
The  probability  of  catching  a  squid  was  significantly 
less  during  maximum  flood/maximum  ebb  (p  <  0.001; 
\|/  =  0.56)  and  late  flood  through  early  ebb  (p  =  0.010; 


\|r  =  0.67)  than  during  late  ebb  through  early  flood 
(Table  1).  The  probability  of  catching  a  squid  increased 
with  depth  to  a  point  after  which  there  was  no 
detectable  difference  in  squid  catch  probability;  there 
was  significantly  higher  probability  of  catching  a  squid 
at  depths  of  5.0  to  9.9  m  than  at  depths  <5.0  m  (p  = 
0.001;  \\t  =  2.16)  and  an  even  higher  probability  of 
catching  a  squid  at  depths  between  10.0  and  14.9  m 
than  at  depths  <5.0  m  (p  <  0.001;  \|/  =  3.62),  but  no  sig¬ 
nificant  difference  in  catch  probability  was  detected 
when  depths  >15.0  m  were  compared  with  depths  <5.0 
m  (p  =  0.316)  (Table  1). 

Salinity  (p  <0.001;  y  =  4.26)  and  water  temperature 
(p  <0.001;  \|/  =  2.42)  had  significant  impacts  on  squid 
catch  probability,  although  the  effects  were  con¬ 
founded  by  a  water  temperature  x  salinity  interaction 
(p  <0.001;  \(/  =  0.91:  Table  1).  Logit  comparisons  of  all 
water  temperature/salinity  combinations  revealed  that 
only  at  low  temperatures  (5  to  9°C)  did  the  probability 
of  catching  at  least  one  squid  not  increase  significantly 
with  increased  salinity  (Fig.  5).  Otherwise,  the  proba¬ 
bility  of  catching  at  least  1  squid  generally  increased 
with  increased  salinity  and  water  temperature.  When 
both  salinity  and  water  temperature  were  high,  there 
was  about  a  70%  probability  of  catching  a  squid  in  a 


Table  1.  LoIIiguncula  brevis.  Multivariate  logistic  regression  table  of  squid  presence/absence  within  the  Chesapeake  Bay.  Five 
discrete  factors  (tide,  depth,  dissolved  oxygen,  year,  and  season),  2  scaled  continuous  factors  (sahnity  and  water  temperature)  and 
1  interaction  (water  temperature  x  salinity)  are  included  in  the  model.  Reference  levels  (r)  and  incremental  units  (u)  are  included 
in  parentheses  next  to  discrete  and  scaled  continuous  factors,  respectively.  The  logistic  regression  model  was  highly  significant: 
-2  log-likehhood  =  -786.15;  goodness-of-fit  statistic,  G  =  665.28,  df  =  15,  p  <0.001.  Goodness-of-fit  (Hosmer-Lemeshow  method)  = 

chi-square  8.804;  df  =  8;  p  =  0.359 


Predicfor 

Coefficient 

Z-value 

p-value 

Odds 
ratio  (V(r) 

Lower  95  % 

Cl 

Upper  95  % 

Cl 

Consfant 

-15.836 

-12.62 

<0.001 

Tide  (r  =  late  ebb  -  early  flood) 

Max.  flood  /  max.  ebb 

-0.577 

-3.63 

<0.001 

0.56 

0.41 

0.77 

Late  flood  -  early  ebb 

-0.402 

-2.59 

0.010 

0.67 

0.49 

0.91 

Depth  (r  =  0.0  -4.9  m) 

5.0  -  9.9  m 

0.771 

3.22 

0.001 

2.16 

1.35 

3.46 

10.0 -14.9  m 

1.286 

5.27 

<0.001 

3.62 

2.24 

5.83 

>15.0  m 

0.302 

1.00 

0.316 

1.35 

0.75 

2.44 

Salinity  (u  =  2  %o) 

1.449 

7.02 

<0.001 

4.26 

2.84 

6.38 

Water  temperature  (u  =  3°C) 

0.884 

6.63 

<0.001 

2.42 

1.86 

3.14 

Dissolved  oxygen  (r  =  <  3  mg  O2  1  *) 

>3  mg  O2  1  ^ 

2.387 

4.49 

<0.001 

10.88 

3.84 

30.88 

Year  (r  =  1996) 

1993 

1.084 

3.86 

<0.001 

2.96 

1.70 

5.13 

1994 

0.927 

3.30 

0.001 

2.53 

1.46 

4.38 

1995 

1.816 

6.87 

<0.001 

6.15 

3.66 

10.32 

1997 

1.266 

4.79 

<0.001 

3.54 

2.11 

5.95 

Season  (r  =  spring) 

Summer 

0.934 

3.57 

<0.001 

2.55 

1.52 

4.25 

Fall 

2.494 

11.07 

<0.001 

12.11 

7.79 

18.83 

Water  temperature  x  Salinity 

-0.099 

-3.49 

<0.001 

0.91 

0.86 

0.96 

probabi'rtV 
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trawl  (Fig.  5).  Based  on  the  odds  ratios  where  catch 
probability  increased  by  a  factor  of  4.26  with  every  2 %o 
rise  (between  17.5  and  37.1  %o)  and  by  a  factor  of  2.42 
with  every  3°C  rise  (between  5  and  30°C),  salinity  had 
a  larger  influence  on  catch  probability  than  water  tem¬ 
perature.  o 

Q- 

Dissolved  oxygen  was  treated  as  a  dichotomous  dis-  * 
Crete  variable  because  the  only  clear  difference  in  ;| 
catch  probability  occurred  between  normoxic  and 
hypoxic  sites  (i.e.,  no  linear  or  quadratic  relationships 
between  the  logit  and  dissolved  oxygen  levels  were 
present).  Not  surprisingly,  catch  probability  increased 
dramatically  (p  <  0.001;  \|r  =  10.88)  when  sampling  in 
normoxic  as  opposed  to  hypoxic  waters. 

Capture  probability  varied  both  annually  and  sea¬ 
sonally.  The  probability  of  catching  a  squid  was  lowest 
in  1996  and  highest  in  1995,  when  the  probability  of 
catching  a  squid  was  6.15  times  greater  than  in  1996 
(Table  1).  There  was  a  2.55  times  greater  probability  of 
catching  a  squid  in  summer  than  in  spring,  and  a  12.11 
times  greater  probability  of  catching  a  squid  in  fall 
than  in  spring  (Table  1). 

There  was  considerable  annual  and  seasonal  variation 
in  physical  conditions  over  the  period  of  this  study.  In 
1995,  when  the  probability  of  catching  a  squid  was  high¬ 
est,  mean  monthly  bottom- water  salinities  for  all  stations 
sampled  did  not  decline  below  18.5%o  (±7.3  SD),  even 


1993  1994  1995  1996  1997 


Month 


- ■ -  Maiiistem -I- Tributaries 

. ♦ .  Mainstem 

—  O —  Tributaries 


Fig.  6.  Mean  (±1  SE)  monthly  bottom-water  salinities  at  sam¬ 
pling  sites  during  VIMS  Trawl  Surveys  from  1993  to  1997. 
Mean  monthly  salinities  for  all  stations  sampled  (Mainstem 
and  Tributaries),  Chesapeake  Bay  mainstem  sampling  sites 
(Mainstem),  and  tributary  sampling  sites  (Tributaries)  are 
plotted;  17.5  %o  line:  lower  physiological  limit  of  brief  squid 
salinity  tolerance;  shaded  areas:  winter  and  summer  months; 
unshaded  areas:  spring  and  fall  months 


Fig.  5.  Lolliguncula  brevis.  Probability  of  catching  at  least  1 
squid  in  a  trawl  plotted  as  a  function  of  salinity  and  water 
temperature.  Generally,  probabilities  are  enhanced  with 
increasing  salinity  and  temperature;  however,  when  water 
temperatures  are  <9°C,  catch  probability  does  not  increase 
with  increasing  salinity 


during  the  whiter  and  spring,  and  mean  monthly  bottom- 
water  salinities  for  mainstem  Chesapeake  Bay  stations 
remained  above  22.9%o  (±3.3  SD)  (Fig.  6).  Furthermore, 
mean  monthly  bottom-water  temperatures  for  aU  stations 
sampled  were  higher  from  June  through  October  of  1995 
than  during  these  months  in  other  years  (Fig.  7) .  In  1996, 
when  the  probability  of  catching  a  squid  was  lowest, 
mean  monthly  bottom-water  salinities  for  all  stations  re¬ 
mained  below  15.3%o  (±  6.4  SD),  and  mean  monthly  bot¬ 
tom-water  salinities  for  bay  mainstem  stations  were 
<20.7%o  (±4.1SD:Fig.  6). 

During  the  winter,  bottom-water  temperatures  were 
routinely  below  8°C,  overall  mean  salinity  often  de¬ 
creased  relative  to  conditions  in  the  fall,  and  dissolved 
oxygen  levels  were  the  highest  of  the  year  (Figs.  6  to  8). 
In  the  spring,  overall  mean  salinities  often  were  below 
17.5  %o,  with  the  notable  exception  of  1995  when  overall 
mean  salinities  did  not  decrease  below  18.5  %o  for  the 
entire  year.  Relative  to  conditions  in  the  winter,  dis¬ 
solved  oxygen  levels  decreased  and  water  tempera¬ 
tures  increased  in  the  spring.  In  the  summer,  salinities 
increased  relative  to  those  in  spring,  water  tempera¬ 
tures  were  clearly  the  highest  of  the  year,  and  dissolved 
oxygen  levels  were  the  lowest  of  the  year.  In  most 
years,  overall  mean  salinity  was  greater  during  the  fall 
than  at  any  other  time  during  the  year.  Relative  to  the 
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Fig.  7.  Mean  (±1  SE)  monthly  bottom-water  temperatures  at 
sampling  sites  during  VIMS  trawl  surveys  from  1993  to  1997. 
Mean  monthly  temperatures  for  all  stations  sampled  (Main- 
stem  and  Tributaries),  Chesapeake  Bay  mainstem  sampling 
sites  (Mainstem),  and  tributary  sampling  sites  (Tributaries) 
are  plotted.  Boxes:  conditions  from  June-October;  shaded 
areas:  winter  and  summer  months;  unshaded  areas:  spring 
and  fall  months.  Since  there  was  little  variability  in  water  tem¬ 
perature  between  sites  during  a  given  month,  standard  error 
bars  are  not  always  visible 

summer,  water  temperatures  decreased  and  dissolved 
oxygen  levels  increased  during  the  fall  (Figs  6  to  8). 

During  the  course  of  this  study,  LoIIiguncula  brevis 
egg  masses  were  occasionally  observed  within  the 
trawl  mesh.  These  masses  were  collected  in  September 
and  October  within  sandy  bottom  habitats  along  the 
eastern  portion  of  the  Chesapeake  Bay  mainstem 
(Fig.  1).  Individual  egg  cases  were  generally  <2.5  cm 
in  length  and  contained  <30  L.  brevis  paralarvae.  Fur¬ 
thermore,  mature  males  and  ripe  females  were 
observed  in  late  July  to  September. 

DISCUSSION 

Salinity,  water  temperature,  and  dissolved  oxygen 
levels  all  play  critical  roles  in  determining  LoIIiguncula 
brevis  distribution  in  Chesapeake  Bay.  The  data  in  this 
study  indicate  that  there  are  salinity  (17.9%o)  and  water 
temperature  (8.1°C)  thresholds  below  which  L.  brevis 
is  absent  from  Chesapeake  Bay.  There  is  also  a  level  of 
dissolved  oxygen  (3  mg  O2  T^)  below  which  L.  brevis 
presence  decreases  dramatically;  however,  L.  brevis  is 
found  in  Chesapeake  Bay  waters  with  as  low  as  1.9  mg 
O2  r^.  Above  salinity  and  water  temperature  thresh- 
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Fig.  8.  Mean  (±1  SE)  monthly  dissolved  oxygen  levels  of  bot¬ 
tom-water  at  sampling  sites  during  VIMS  trawl  surveys  from 
1993  to  1997.  Mean  monthly  dissolved  oxygen  levels  for  aU  sta¬ 
tions  sampled  (Mainstem  and  Tributaries),  Chesapeake  Bay 
mainstem  sampling  sites  (Mainstem),  and  tributary  sampling 
sites  (Tributaries)  are  plotted.  Shaded  areas:  winter  and  sum¬ 
mer  months;  unshaded  areas:  spring  and  fall  months 

olds,  the  probability  of  catching  a  squid  increases 
steadily  with  increased  salinity  and  water  temperature 
up  to  conditions  of  35  %o  and  29.6°  C,  respectively.  In 
normoxic  waters  (3.0  to  14.6  mg  O2  1”^),  as  opposed  to 
hypoxic  waters  (<3.0  mg  O2  1“*),  the  probability  of 
catching  a  squid  increases  by  a  factor  of  10.  Conse¬ 
quently,  seasonal  and  annual  variation  in  one  or  more 
of  these  factors  has  a  significant  impact  on  squid  distri¬ 
bution  in  Chesapeake  Bay. 

The  physical  limits  detected  in  this  study  are  consis¬ 
tent  with  those  in  previous  studies  on  LoIIiguncula  bre¬ 
vis,  and  provide  further  evidence  that  this  euryhaline 
squid  has  exceptional  tolerance  limits  relative  to  other 
cephalopods.  In  controlled  laboratory  studies,  Hendrix 
et  al.  (1981)  determined  that  L.  brevis  tolerates  salini¬ 
ties  down  to  17.5%o,  a  limit  close  to  that  detected  in  the 
present  study  and  below  which  L.  brevis  demonstrates 
signs  of  hypoosmotic  stress.  The  respiratory  pigment 
of  L.  brevis,  hemocyanin,  is  insensitive  to  salinity 
changes,  has  high  oxygen  affinity,  and  has  little  pH 
dependence  over  a  broad  range  of  salinities  —  charac¬ 
teristics  that  are  important  for  life  in  euryhaline  waters. 
Conversely,  hemocyanin  in  stenohaline  squid,  e.g. 
Loligo  pealei,  is  sensitive  to  salinity  changes,  has  rela¬ 
tively  low  oxygen  affinity,  and  has  large  pH  depen¬ 
dence,  preventing  these  squids  from  entering  low- 
salinity  areas  for  extended  periods  (Mangum  1991). 
The  salinity  limit  determined  in  the  present  study  is  in 
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reasonable  agreement  with  the  field-salinity  limits  of 
18.2 %o  determined  by  Dragovich  &  Kelly  (1967)  along 
the  western  coast  of  Florida  and  20.0 %o  determined  by 
Hixon  (1980)  in  the  northwest  Gulf  of  Mexico.  A  lower 
limit  of  8.5  %o  was  reported  by  Laughlin  &  Livingston 
(1982)  in  an  estuary  in  north  Florida,  however.  The 
low-end  salinity  thresholds  reported  in  the  present 
study  and  in  previous  studies  are  lower  than  those  doc¬ 
umented  for  other  squids  in  the  north  Atlantic,  such  as 
L.  pealei,  L.  plei,  and  Illex  illecebrosus,  which  are  typ¬ 
ically  found  in  waters  >30%o  (O'Dor  et  al.  1977,  Flixon 
1980). 

In  the  Gulf  of  Mexico,  Dragovich  &  Kelly  (1967), 
Hixon  (1980)  and  Laughlin  St  Livingston  (1982)  cap¬ 
tured  Lolliguncula  brevis  in  waters  of  12.6  to  31.6°C, 
11  to  31° C,  and  15  to  31° C,  respectively.  The  results  of 
the  present  study  reveal  that  L.  brevis  may  be  found  in 
slightly  cooler  waters  near  the  northern  limit  of  its 
range.  In  a  study  on  the  respiratory  and  cardiac  func¬ 
tion  of  L.  brevis  performed  in  circulating  and  non-cir¬ 
culating  respirometers.  Wells  et  al.  (1988)  determined 
that  L.  brevis  functions  well  at  water  temperatures  14 
to  27° C.  Although  Wells  et  al.  did  not  report  on  the 
physiology  of  L.  brevis  below  14°  C,  they  did  find  that 
L.  brevis  demonstrates  clear  signs  of  stress,  such  as 
decreased  heartbeat  frequency  and  increased  stroke 
volume,  at  higher  water  temperatures  (27  to  30°  C)  — 
temperatures  at  which  the  squid  were  frequently 
caught  during  the  present  study.  Wells  et  al.  attributed 
the  ability  of  L.  brevis  to  tolerate  water  temperatures  up 
to  31°C  in  nature  to  long  acclimation  periods.  Loligo 
plei  has  a  wide  range  of  temperature  tolerances  (13  to 
30°C),  hke  that  of  Lolliguncula  brevis,  but  most  squids 
have  more  restricted  temperature  ranges,  e.g.  Illex iUe- 
cebrosus  and  Loligo  pealei,  which  are  found  predomi¬ 
nantly  in  waters  of  8  to  15°C  and  7  to  22°C,  respectively 
(Hixon  1980,  Whitaker  1980,  Roper  et  al.  1984). 

The  discovery  of  Lolliguncula  brevis  in  hypoxic 
waters  of  Chesapeake  Bay  is  consistent  with  previous 
studies,  and  distinguishes  L.  brevis  from  other  squids. 
Vecchione  &  Roper  (1991)  observed  L.  brevis  in 
hypoxic  waters  (0.7  mg  O2  1“^)  of  the  western  North 
Atlantic  using  remotely  operated  submersibles,  and 
Vecchione  (1991b)  documented  L.  brevis  in  hypoxic 
waters  (<3  mg  O2  1“*)  within  coastal  and  estuarine 
waters  of  southwestern  Louisiana  during  trawl  sur¬ 
veys.  The  unique  characteristics  of  the  hemocyanin  of 
L.  brevis,  especially  its  high  O2  affinity,  and  the  ability 
of  L.  brevis  to  build  up  large  oxygen  debts  (unlike 
other  coastal  squids)  are  beneficial  for  L.  brevis  in 
waters  with  low  dissolved  oxygen  levels  (Mangum 
1991,  Finke  et  al.  1996).  Wells  et  al.  (1988)  determined 
that  L.  brevis  activity  levels  generally  decrease  when 
dissolved  oxygen  levels  are  <3  mg  O2  1“^,  but  they 
found  that  some  L.  brevis  show  no  behavioral  changes 


down  to  levels  of  2  mg  O2  1“^.  Because  energy  avail¬ 
ability  for  cellular  ATPases  declines  rapidly  under 
hypoxic  conditions,  Finke  et  al.  (1996)  and  Zielinski  et 
al.  (2000)  suggested  that  L.  brevis  may  spend  only  brief 
periods  in  hypoxic  waters.  Excursions  by  L.  brevis  into 
hypoxic  waters,  however  short  in  duration,  may  be 
performed  to  exploit  a  food  niche  or  avoid  predation  as 
proposed  by  Vecchione  (1991a).  Two  possible  fishes 
upon  which  L.  brevis  may  prey  in  hypoxic  waters  (pro¬ 
vided  they  are  within  the  proper  size  range  [i.e.,  length 
of  squid  arms  or  less])  are  the  spot  Leiostomus  xanthu- 
rus  and  the  Atlantic  croaker  Micropongonias  undula- 
tus — the  fishes  most  commonly  encountered  in  such 
waters  during  VIMS  Trawl  Surveys 

The  synergistic  and  independent  effects  of  salinity, 
water  temperature  and  dissolved  oxygen,  which  all 
increase  catch  probabilities  when  elevated,  drive  sea¬ 
sonal  and  annual  distribution  patterns  in  Chesapeake 
Bay.  Although  high  dissolved  oxygen  levels  were 
observed  throughout  Chesapeake  Bay  during  the  win¬ 
ter,  salinity  decreased  (relative  to  conditions  in  the  fall) 
and  water  temperature  dropped  below  8°C,  triggering 
an  exodus  of  Lolliguncula  brevis  from  the  bay.  When 
water  temperatures  increased  to  approximately  10°C 
in  April  and  May  (spring),  squid  re-entered  Chesa¬ 
peake  Bay  in  areas  where  salinity  was  >17.5%o.  Similar 
migratory  responses  of  L.  brevis  to  water  temperatures 
were  reported  by  Hixon  (1980),  who  found  that  L.  bre¬ 
vis  move  offshore  from  Galveston  Bay,  Texas,  in 
December  through  February,  when  mean  water  tem¬ 
peratures  are  12.8  to  13.9°C,  and  re-enter  inshore  bay 
habitats  in  March,  when  temperatures  are  above  15°C. 
The  extent  to  which  L.  brevis  enters  Chesapeake  Bay 
in  the  spring  and  throughout  the  year  depends  heavily 
on  salinity.  In  years  when  salinity  was  low  throughout 
the  year,  e.g.  1996,  squid  entered  Chesapeake  Bay 
later  (May  as  opposed  to  April),  were  captured  less  in 
trawls,  and  were  found  exclusively  in  the  lower,  high- 
salinity  regions  of  the  bay.  In  years  when  salinity  was 
high  in  the  spring  and  throughout  the  year,  e.g.  1995, 
squid  entered  Chesapeake  Bay  early  (April  as  opposed 
to  May),  were  more  prevalent  in  trawls,  and  were  cap¬ 
tured  in  a  wider  region  of  the  bay.  High  water  temper¬ 
atures  from  June  through  October  also  contributed  to 
the  observed  high  catch  probabilities  in  1995. 

High  salinity  and  water  temperatures  or  high  salinity 
alone  contributed  to  greater  squid  catch  probabilities 
in  the  summer  and  fall  than  in  the  winter  and  spring. 
Given  that  squid  catch  probabilities  were  highest  in 
the  fall,  when  waters  were  less  hypoxic,  often  more 
saline,  but  cooler  relative  to  waters  during  the  summer, 
it  appears  that  the  benefits  of  high  salinity  and  more 
widespread  normoxia  outweigh  the  disadvantages  of 
lower  water  temperatures.  This  is  not  surprising  con¬ 
sidering  the  fact  that  salinity  had  a  greater  impact  on 
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squid  catch  probability  than  water  temperature  in  the 
logistic  regression  model,  and  hypoxic  water  lowers 
catch  probabilities  substantially.  Laughlin  &  Livingston 
(1982)  also  found  that  the  spatial  distribution  of  brief 
squid  within  a  Florida  estuary  is  strongly  influenced  by 
salinity  as  well  as  habitat  structure,  while  temperature 
plays  a  lesser  role.  Interestingly,  in  non-estuary  set¬ 
tings  where  salinity  does  not  fluctuate  as  significantly, 
there  is  considerably  less  seasonal  variation  in  the  fre¬ 
quency  of  occurrence  of  loliginids  (Whitaker  1980). 

The  regions  where  squid  were  encountered  most 
frequently  (i.e.,  central  Chesapeake  Bay  near  the 
mouth  and  along  eastern  portions  of  the  bay  where 
waters  are  deepest  and  heavily  influenced  by  intrud¬ 
ing  salt  wedges)  typically  have  the  highest  salinities 
throughout  the  year.  Consequently,  salinity  assuredly 
played  an  important  role  in  the  observed  distribution 
patterns;  however,  depth  played  a  significant  role  as 
well.  Brief  squid  are  shallow-water  cephalopods,  pre¬ 
ferring  habitats  <30  m  in  Galveston  Bay,  Texas  (LaRoe 
1967,  Hixon  1980)  and  along  the  eastern  US  coast,  as 
shown  in  the  NMFS/NEFSC  surveys.  The  bathymetry 
of  Chesapeake  Bay  is  predominantly  <  20  m,  and  con¬ 
sequently  the  bay  provides  optimal  depths  for  squid 
residence.  The  dafa  presented  is  this  study  indicate 
that  there  are  further  preferences  at  depths  less  than 
30  m  within  the  Chesapeake  Bay  ecosystem,  with 
squid  favoring  depths  of  10  to  15  m  even  over  deeper 
waters,  which  are  more  saline  but  also  more  hypoxic. 
These  depths  are  found  in  the  central  channels  of 
Chesapeake  Bay — areas  where  the  highest  overall 
number  of  squid  were  encountered.  Therefore,  when 
salinities  and  water  temperatures  are  within  tolerance 
limits,  brief  squid  prefer  central  channel  areas  10  to  15 
m  in  depth.  Laughlin  &  Livingston  (1982)  found  a  simi¬ 
lar  preference  for  channels  by  brief  squid  in 
Appalachicola  estuary,  Florida. 

The  higher  catch  probability  observed  during  late 
ebb  through  early  flood  compared  with  maximum 
flood/maximum  ebb  or  late  flood  through  early  ebb 
raises  some  interesting  questions.  This  result  may  be, 
in  part,  a  product  of  a  reduction  in  water  volume  within 
the  Chesapeake  Bay  ecosystem  and  the  consequent 
concentration  of  squid  communities,  resulting  in  easier 
entrapment  of  squid  in  trawls.  Llowever,  it  also  proba¬ 
bly  reflects  a  behavioral  preference.  LoIIiguncula  bre¬ 
vis  swims  at  low  speeds;  it  has  a  critical  swimming 
speed  of  <25  cm  s“^,  and  demonstrates  parabolic  oxy¬ 
gen  consumption  curves  as  a  function  of  speed,  wifh 
oxygen  consumption  minima  at  low/intermediate 
speeds  (6  to  12  cm  s“^:  Bartol  et  al.  2001  a,b).  Conse¬ 
quently,  lower  currents  associated  with  late  ebb 
though  early  flood  tides  may  be  advantageous  for 
maneuvering  in  the  bay  environment  and  for  keeping 
pace  with  the  oncoming  current  (hovering),  which  is 


frequently  observed  by  squids  in  nature  (O'Dor  et  al. 
1994).  Given  that  Laughlin  &  Livingston  (1982)  found  L. 
brevis  to  be  more  abundant  in  deep  channels  where 
currents  are  higher  than  in  shallow  (<4  m),  enclosed 
regions,  there  may  be  a  lower  limit  to  flow  velocity,  with 
L.  brevis  seeking  out  currents  matching  optimal  swim¬ 
ming  speeds,  as  reported  for  other  squids  (O’Dor  et  al. 
2001).  One  important  question  not  addressed  directly  in 
the  study  is:  where  do  squid  go  during  tidal  cycles  when 
catch  probabilities  are  low?  Given  that  trawls  conducted 
in  this  survey  sampled  only  the  lower  water  column  of 
fhe  Chesapeake  Bay  system,  it  is  conceivable  that  squid 
swam  higher  in  the  water  column  during  the  different 
tidal  phases  or  simply  swam  out  to  sea.  Future  study, 
possibly  involving  the  deployment  of  autonomous 
underwater  vehicles  (AUVs),  is  necessary  to  determine 
squid  movements  and  behaviors  within  Chesapeake 
Bay  as  a  function  of  tidal  phase  and  flow  speed. 

Brief  squid  are  clearly  abundant  in  Chesapeake  Bay 
relative  to  other  nekton  captured  in  the  VIMS  trawl 
surveys  from  1993  to  1997,  ranking  in  the  upper  6  to 
12%  of  nekfonic  catches  in  terms  of  total  number  of 
organisms  capfured  in  3  of  the  5  sampling  years  and  in 
the  upper  20%  in  all  years  but  1996.  Moreover,  during 
the  fall  of  some  years,  squid  ranked  second  to  bay 
anchovies.  The  results  of  this  study  indicate  that  2  size 
classes  of  LoIIiguncula  brevis  use  Chesapeake  Bay: 
one  class  (<60  mm  DML)  uses  the  bay  in  the  spring 
through  the  fall,  and  a  second,  less  abundant  class 
(>70  mm  DML)  is  most  prevalent  in  the  spring.  The 
absence  of  a  distinct  intermediate  size  class  (~60  mm 
DML)  in  Chesapeake  Bay  in  the  fall  (or  at  any  other 
time  during  the  year),  as  was  detected  in  the 
NMFS/NEFSC  surveys  performed  along  the  Atlantic 
coast,  suggests  that  L.  brevis  may  leave  the  bay  once  a 
specific  size  is  reached  and  head  for  coastal  waters  of 
the  Atlantic  Ocean.  Fleavy  utilization  of  cerfain  habi- 
fats  by  smaller  individuals  and  subsequent  migration 
of  larger  individuals  is  common  among  cephalopods 
(Whitaker  1980,  Hixon  1983,  Nesis  1983,  Okutani 
1983,  Summers  1983,  Nagasawa  et  al.  1993,  Hanlon  & 
Messenger  1996). 

Limited  field  observations  suggest  that  larger  brief 
squid  may  be  using  Chesapeake  Bay  as  spawning 
grounds,  at  least  during  late  July  through  early  Sep¬ 
tember.  This  is  consistent  with  the  study  of  Vecchione 
(1982),  who  found  paralarvae  near  the  mouth  of 
Chesapeake  Bay  during  the  warmest  months  of  late 
summer.  Within  the  Gulf  of  Mexico,  spawning  occurs 
to  some  degree  year-round,  with  peak  spawning 
occurring  during  April  to  July  and  September  to 
November  (LaRoe  1967,  Hixon  1980,  Jackson  et  al. 
1997).  More  work  on  fecundity  and  paralarval  ecology 
needs  to  be  performed  in  Chesapeake  Bay  to  deter¬ 
mine  if  spawning  occurs  in  the  spring  (when  larger 
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Lolliguncula  brevis  are  most  abundant),  summer,  and 
fall,  and  if  fhere  are  peak  spawning  evenfs  throughout 
the  year.  Given  that  statolith  studies  have  revealed 
that  L.  brevis  live  approximately  100  to  200  d  in  the 
Gulf  of  Mexico  (Jackson  ef  al.  1997),  it  is  likely  that 
spawning  occurs  to  some  degree  most  of  fhe  fime  that 
brief  squid  are  in  Chesapeake  Bay.  Early  juveniles  that 
either  hatch  in  Chesapeake  Bay  or  that  migrate  into 
the  bay  from  offshore  habifats  may  remain  in  the  bay 
to  exploit  its  rich  food  supply,  i.e.,  crustaceans  and 
juvenile  finfishes. 

The  common  occurrence  of  brief  squid  within  Chesa¬ 
peake  Bay,  a  euryhaline,  highly  variable  environment 
that  is  home  to  hundreds  of  species  of  fishes  each  year, 
is  intriguing.  Squid  in  general  are  thought  to  be  ham¬ 
pered  by  a  deficient  physiology  and  are  notorious  for 
possessing  an  inefficienf  means  of  propulsion.  Conse- 
quenfly,  squid  are  thought  to  have  evolved  in  environ¬ 
ments  such  as  deep-sea  or  offshore,  pelagic  habitafs 
where  inferaction  with  more  efficient  nekton  is  mini¬ 
mized  (O'Dor  &  Webber  1986,  Wells  1994).  Nonethe¬ 
less,  brief  squid  are  present  in  an  inshore,  physiologi¬ 
cally  taxing  habitat  rich  with  nektonic  fauna.  Furfher 
research  info  fhe  physiology  and  locomotive  abilities  of 
brief  squid,  especially  the  smaller  individuals  that 
appear  to  be  most  abundant  within  Chesapeake  Bay, 
may  further  our  evolutionary  and  ecological  under¬ 
standing  of  this  unique  cephalopod. 
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